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ABSTRACT 

The observed 21-cm signal from the epoch of reionization will be distorted along 
the line-of-sight by the peculiar velocities of matter particles. These redshift-space 
distortions will affect the contrast in the signal and will also make it anisotropic. 
This anisotropy contains information about the cross-correlation between the matter 
density field and the neutral hydrogen field, and could thus potentially be used to 
extract information about the sources of reionization. In this paper, we study a collec¬ 
tion of simulated reionization scenarios assuming different models for the sources of 
reionization. We show that the 21-cm anisotropy is best measured by the quadrupole 
moment of the power spectrum. We find that, unless the properties of the reionization 
sources are extreme in some way, the quadrupole moment evolves very predictably 
as a function of global neutral fraction. This predictability implies that redshift-space 
distortions are not a very sensitive tool for distinguishing between reionization sources. 
However, the quadrupole moment can be used as a model-independent probe for con¬ 
straining the reionization history. We show that such measurements can be done to 
some extent by first-generation instruments such as LOFAR, while the SKA should be 
able to measure the reionization history using the quadrupole moment of the power 
spectrum to great accuracy. 

Key words: cosmology:dark ages, reionization, first stars—methods: numerical 


1 INTRODUCTION 

One of the periods in the history of our Universe 
about which we know the least is the epoch of reion¬ 
ization (EoR). During this epoch the first sources of 
light formed and gradually ionized the neutral hydrogen 
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(Hi) in the intergalactic medium (IGM). Our present un¬ 
derstanding of this epoch is mainly constrained by ob¬ 
servations of the cosmic microwave background radia¬ 
tion (CMBR) (Komatsu et al. 2011; Planck Collaboration 
2015) and the absorption spectra of high redshift quasars 
(Becker et al. 2001; Fan et al. 2003; White et al. 2003; 
Goto et al. 2011; Becker et al. 2015). These observations 
suggest that reionization was an extended process, span¬ 
ning over the redshift range 6 < 2 : < 15 (see 
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e.g. Alvarez et al. 2006; Mitra, Ferrara & Choudhury 2013; 
Mitra, Choudhury & Ferrara 2015; Robertson et al. 2015; 
Bouwens et al. 2015). Such indirect observations are, how¬ 
ever, limited in their ability to answer several important 
questions regarding the EoR. These unresolved issues in¬ 
clude the precise duration and timing of reionization, the 
properties of major ionizing sources, the relative contribu¬ 
tion to the ionizing photon budget from various kinds of 
sources, and the typical size and distribution of ionized bub¬ 
bles. 

Observations of the redshifted 21-cm line, originating 
from spin flip transitions in neutral hydrogen atoms, are ex¬ 
pected to be the key to resolving many of these long standing 
issues. The brightness temperature of the redshifted 21-cm 
line directly probes the H l distribution at the epoch where 
the radiation originated. Observing this line enables us, in 
principle, to track the entire reionization history as it pro¬ 
ceeds with redshift. 

Motivated by this, a huge effort is underway to 
detect the redshifted 21-cm signal from the EoR us¬ 
ing low frequency radio interferometers, such as the 
GMRT (Paciga et al. 2013), LOEAR (Yatawatta et al. 
2013; van Haarlem et al. 2013; Jelic et al. 2014), MWA 
(Tingay et al. 2013; Bowman et al. 2013), PAPER 
(Parsons et al. 2014) and 21CMA (Wang et al. 2013). 
These observations are complicated to a large degree by 
foreground emissions, which can be ~ 4 — 5 orders of magni¬ 
tude stronger than the expected signal (e.g. Di Matteo et al. 
2002; Ali, Bharadwaj & Chengalur 2008; Jelic et al. 2008), 
and system noise (Morales 2005; McQuinn et al. 2006). 
So far only weak upper limits on the 21-cm signal have 
been obtained (Paciga et al. 2013; Dillon et al. 2014; 
Parsons et al. 2014; Ali et al. 2015). 

Owing to the low sensitivity of the first generation 
interferometers, they will probably not be capable of di¬ 
rectly imaging the Hi distribution. This will have to wait 
for the arrival of the extremely sensitive next genera¬ 
tion of telescopes such as the SKA (Mellema et al. 2013, 
2015; Koopmans et al. 2015). The first generation telescopes 
are instead expected to detect and characterize the sig¬ 
nal through statistical estimators such as the variance (e.g. 
Patil et al. 2014) and the power spectrum (e.g. Pober et al. 
2014). 

Many studies to date have focused on the spherically 
averaged 21-cm power spectrum (e.g. McQuinn et al. 2007; 
Lidz et al. 2008; Barkana 2009; Iliev et al. 2012). By aver¬ 
aging in spherical shells in Fourier space, one can obtain 
good signal-to-noise for the power spectrum, while still pre¬ 
serving many important features of the signal. However, the 
spherically averaged power spectrum does not contain all 
information about the underlying 21-cm field. For example, 
as reionization progresses, the fluctuations in the Hi field 
will make the 21-cm signal highly non-Gaussian. This non- 
Gaussianity can not be captured by the power spectrum, 
but requires higher-order statistics (Bharadwaj & Pandey 
2005; Mellema et al. 2006; Watkinson & Pritchard 2014; 
Mondal, Bharadwaj & Majumdar 2015). 

Another important effect that can not be captured 
fully by the spherically averaged power spectrum is 
the effect of redshift-space distortions caused by the 
peculiar velocities of matter. The coherent inflows of 
matter into overdense regions and the outflows of mat¬ 


ter from underdense regions will produce an additional 
red- or blueshift in the 21-cm signal on top of the cos¬ 
mological redshift, changing the contrast of the 21-cm 
signal, and making it anisotropic (Bharadwaj & Ali 
2004; Barkana & Loeb 2005; McQuinn et al. 2006; 
Mao et al. 2012; Majumdar, Bharadwaj & Choudhury 
2013; Jensen et al. 2013). The redshift-space distortions in 
the 21-cm signal have previously generated interest due to 
the possibility of extracting the purely cosmological matter 
power spectrum (Barkana & Loeb 2005; McQuinn et al. 
2006; Shapiro et al. 2013). However, they also carry 
interesting astrophysical information. 

It has been shown in previous studies 
(Majumdar, Bharadwaj & Choudhury 2013; Jensen et al. 
2013; Ghara, Choudhury & Datta 2015) that the anisotropy 
in the 21-cm signal due to the redshift-space distortions 
will depend on the topology of reionization, or in other 
words on the properties of the sources of reionization. In 
this paper, we explore the prospects of using the 21-cm 
redshift-space distortions to obtain information about the 
sources of reionization. We calculate the evolution of the 
power spectrum anisotropy for a collection of simulated 
reionization scenarios assuming different properties of the 
sources of ionizing photons. We also investigate different 
ways of quantifying this anisotropy, and study the prospects 
of distinguishing between the source models in observations 
with current and upcoming interferometers. 

The structure of this paper is as follows. In Section 

2, we describe the simulations and source models that we 
have used to generate a collection of reionization scenarios. 
Here, we also describe the method that we use to include the 
redshift-space distortions in the simulated signal. In Section 

3, we evaluate different methods of quantifying the power 
spectrum anisotropy and discuss the interpretation of the 
anisotropy using the quasi-linear model of Mao et al. (2012). 
In Section 4, we describe the observed large scale features 
of the redshift-space anisotropy for the different reionization 
scenarios that we have considered. In Section 5, we discuss 
the feasibility of observing these anisotropic signatures in 
present and future radio interferometric surveys, and what 
we can learn from the redshift- space distortions about the 
reionization history. Finally, in Section 6 we summarise our 
findings. 

Throughout the paper we present our results for the 
cosmological parameters from WMAP five year data re¬ 
lease h = 0.7, fim = 0.27, Qa = 0.73, = 0.0226 

(Komatsu et al. 2009). 


2 REIONIZATION SIMULATIONS FOR 
DIFFERENT SOURCE MODELS 

To study the effects of the sources of reionization on the 
redshift-space anisotropy of the 21-cm signal, we simulate 
a number of different reionization scenarios, assuming dif¬ 
ferent source properties. All of these simulations are based 
on a single A-body simulation of the evolving dark matter 
density field. 
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2.1 A'^-body simulations 

The A^-body simulations were carried out as part 
of the PRACE4LOFAR project (PRACE projects 
2012061089 and 2014102339) with the cubep^m code 
(Harnois-Deraps et al. 2013), which is based on the older 
code PMFAST (Merz, Pen & Trac 2005). Gravitational forces 
are calculated on a particle-particle basis at close distances, 
and on a mesh for longer distances. The size of our simula¬ 
tion volume was 500/h = 714 Mpc (comoving) along each 
side. We used 6912^ particles of mass 4.0 x 10^ Mq on a 
13824® mesh, which was then down-sampled to a 600® grid 
for modelling the reionization. Further details of the A-body 
simulation can be found in Dixon et al. (2015) (in prepara¬ 
tion). 

For each redshift output of the A-body simulation, 
haloes were identified using a spherical overdensity scheme. 
The minimum halo mass that we have used for our reioniza¬ 
tion simulations is 2.02 x 10® Mq. 


2.2 Reionization simulations 

To generate the ionization maps for our different reion¬ 
ization scenarios, we used a modified version of the semi- 
numerical code described in Choudhury, Haehnelt & Regan 
(2009); Majumdar, Bharadwaj & Choudhury (2013); 
Majumdar et al. (2014). This method has been tested 
extensively against a radiative transfer simulation and it 
is expected to generate the 21-cm signal from the EoR 
with an accuracy of > 90% (Majumdar et al. 2014) for the 
length scales that we will deal with in this paper. 

Like most other semi-numerical methods for 
simulating the EoR 21-cm signal, our scheme is 
based on the excursion-set formalism developed by 
Furlanetto, Zaldarriaga & Hernquist (2004), making it 
similar to the methods described by Zahn et al. (2007); 
Mesinger & Furlanetto (2007); Santos et al. (2010). In this 
method, one compares the average number of photons 
in a specific volume with the average number of neutral 
hydrogen atoms in that volume. Here we assume that the 
neutral hydrogen follows the dark matter distribution. 

Once we generate the dark matter density Held at a 
fixed redshift using the A-body simulation, then depending 
on the source model of the reionization scenario under con¬ 
sideration (that we will discuss later in this section), we also 
generate an instantaneous ionizing photon field at the same 
redshift in a grid of the same size as that of the dark matter 
density field. These density fields of photons and neutral hy¬ 
drogen are generally constructed in a grid coarser than the 
actual A-body resolution. In our case the resolution of this 
grid is 1.19 Mpc, i.e. 600® cells. 

Next, to determine the ionization state of a grid cell, 
we compare the average number density of ionizing photons 
and neutral hydrogen atoms around it within a spherical 
volume. The radius of this smoothing sphere is then gradu¬ 
ally increased, starting from the grid cell size and going up 
to the assumed mean free path® of the photons at that red¬ 
shift. If, for any radius of this smoothing sphere, the average 

® The mean free path of ionizing photons at high redshifts is 
largely an unknown quantity till date. We have used a fixed max¬ 
imum smoothing radius of 70 comoving Mpc at all redshifts for 


number density of photons becomes greater than or equal 
to the average number density of the neutral hydrogen, the 
cell is flagged as ionized. The same procedure is repeated for 
all the grid cells and an ionization field at that specific red¬ 
shift is generated. A more detailed description of this simula¬ 
tion method can be found in Choudhury, Haehnelt & Regan 
(2009) and Majumdar et al. (2014). 


2.2.1 Source Models 

Each of our simulated reionization scenarios consists of a 
different combination of various sources of ionizing photons. 
Here, we describe these different source types, and in Section 

2.2.2 we describe the reionization scenarios. 


(i) Ultraviolet photons from galaxies (UV): 

In most reionization models, galaxies residing in the col¬ 
lapsed dark matter halos are assumed to be the major 
sources of ionizing photons. To date, much is unknown about 
these high-redshift galaxies and the characteristics of their 
radiation. Thus most simulations assume that the total num¬ 
ber of ionizing photons contributed by a halo of mass Mj, 
which is hosting such galaxies is simply: 


A.,(Mh) = N, 


M\Slh 

mpllni 


( 1 ) 


Here, Aion is a dimensionless constant which effectively rep¬ 
resents the number of photons entering in the IGM per 
baryon in collapsed objects and rup is the mass of a pro¬ 
ton or hydrogen atom. Reionization simulations (whether 
radiative transfer or semi-numerical) which adopt this kind 
of model for the production of the major portions of 
their ionizing photons generally produce a global “inside- 
out” reionization scenario (see e.g. Mellema et al. 2006; 
Zahn et al. 2007; Ghoudhury, Haehnelt & Regan 2009; 
Mesinger, Furlanetto & Gen 2011 etc.). 

We assume that all ionizing photons generated by these 
kinds of sources® are in the ultraviolet part of the spectrum. 
Thus, they only affect the IGM locally, up to a distance 
limited by their mean free path. 

(ii) Uniform ionizing background (UIB): 

The observed population of galaxies at high redshifts, 
seems unable to keep the universe ionized, unless there 
is a significant increase in the escape fraction of the 


our fiducial reionization sources (i.e. UV sources), which is con¬ 
sistent with the findings of Songaila & Cowie (2010) at z ~ 6. 

^ Note that in our reionization scenarios we allow only halos 
of mass > 2.02 X 10® Mq to host UV photon sources. In re¬ 
ality faint galaxies hosted by low mass halos (10® < < 

10® Mq) may produce significant amounts of UV photons dur¬ 
ing the EoR. However, there is a possibility that star forma¬ 
tion in these low mass halos may also get suppressed once they 
have been ionized and heated to temperatures of ~ 10^ K (see 
e.g. Couchman & Rees 1986; Gnedin 2000; Dijkstra et al. 2004; 
Okamoto, Gao & Theuns 2008). This poses an uncertainty in 
their actual role during the EoR. Further, even when one includes 
them in the simulations, it is very unlikely that the resulting ion¬ 
ization and 21-cm topology will depart from its global “inside-out” 
nature, and the anisotropy of the 21-cm power spectrum—which 
is the focus of this paper—will not be severely sensitive to the 
presence of low mass sources (we discuss this in further details in 
Appendix A). 
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ionizing photons from them with the increasing redshift 
or the galaxies below the detection thresholds of the 
present day surveys contribute a significant fraction of 
the total ionizing photons (Kuhlen & Faucher-Giguere 2012; 
Mitra, Ferrara & Choudhury 2013). An alternative possibil¬ 
ity is that, if sources of hard X-ray photons (such as active 
galactic nuclei or X-ray binaries) were common in the early 
Universe, these could give rise to a more or less uniform 
ionizing background. Hard X-rays would easily escape their 
host galaxies and travel long distances before ionizing hy¬ 
drogen (McQuinn 2012; Mesinger, Ferrara & Spiegel 2013). 

We model this type of source as a completely uniform 
ionizing background that provides the same number of ion¬ 
izing photons at every location. The extreme case of a 100% 
contribution of the ionizing photons from this kind of back¬ 
ground would lead to a global “outside-in” reionization. 

(iii) Soft X-ray photons (SXR): 

The escape fraction of extreme ultraviolet photons from 
their host galaxies is a hotly debated issue. Soft X-ray pho¬ 
tons, on the other hand, would have little difficulty escap¬ 
ing from their host sources into the IGM. Also, there is a 
possibility that the X-ray production could have been more 
prevalent in the high redshift galaxies than their low redshift 
counterparts (Mirabel et al. 2011; Fragos et al. 2013). 

When including this type of photons in our simulations 
we have considered them to be uniformly distributed around 
their source halos within a radius equal to their mean free 
path at that redshift. To estimate the mean free path of these 
soft X-ray photons we have used equation (1) in McQuinn 
(2012) which is dependent on the redshift of their origin 
and the frequency of the photon. For simplicity, we assume 
that all the halos that we have identified as sources will 
produce soft X-ray photons and all of these photons will 
have the same energy (200 eV). A significant contribution by 
these soft X-ray photons may lead to a more homogeneous 
reionzation scenario than the case when the major portion 
of photons are in the ultraviolet. 

(iv) Power law mass dependent efficiency (PL): 

For source type (i) we have assumed that the number of 

UV photons generated by a galaxy is proportional to its host 
halo mass. Instead of this it can also be assumed that the 
number of photons contributed by a halo follows a power 
law iV^(Mh) oc M^. 

We consider two cases, where the power law index is 
equal to 2 and 3. With this source model, higher-mass 
haloes produce relatively more ionizing photons, giving rise 
to fewer, but larger ionized regions. This model provides 
a crude approximation of a situation in which reionization 
is powered by rare, bright sources, such as quasars. While 
the number density of quasars at high redshifts is not ex¬ 
pected to be high enough to drive reionization (see e.g. 
Madau, Haardt & Rees 1999), this model serves as an ex¬ 
treme illustration of the type of topology expected from rare 
ionizing sources. 

2.2.2 Reionization scenarios 

The aim of this paper is to study the effect of the different 
source types described above on the 21-cm power spectrum 
anisotropy. To do this, we construct several reionization sce¬ 
narios by combining the source types in different ways, as 
summarized in Table 1. 



Figure 1. The evolution of the mass averaged neutral fraction 
with redshift for our fiducial reionization scenario. We tune all 
other reionization scenarios to follow the same reionization his¬ 
tory. 


In our fiducial reionization scenario, 100% of the ion¬ 
izing photons come from galaxies residing in dark matter 
halos of mass > 2.02 x 10® Mq , resulting in a global Inside- 
out reionization topology. The evolution of the reionization 
source population in this scenario thus follows the evolu¬ 
tion of the collapsed fraction. Reionization starts around 
2 ; ~ 16 (when the first ionization sources formed) and hy z = 
7.221 the mass averaged neutral fraction (designated by xhi 
throughout this paper) of the IGM reaches ~ 0.2, as shown 
in Figure 1. At this point, the 21-cm signal is likely too 
weak to detect, at lerst with first-generation interferometers 
(Jensen et al. 2013; Majumdar, Bharadwaj & Ghoudhury 
2013; Datta et al. 2014; Majumdar et al. 2014). 

We tune Won at each redshift for all our reionization 
scenarios to follow the same evolution of xhi with z as the 
fiducial model (see Figure 1). This ensures that any differ¬ 
ences in the 21-cm signal and its anisotropy across different 
reionization scenarios are due solely to the source types, and 
not due to the underlying matter distribution. Of course this 
will force all scenarios except the fiducial one to follow an ar¬ 
tificial reionization history. However, our aim here is not to 
produce the most realistic reionization scenarios, but rather 
to produce a wide range of reionization topologies and thus 
a significant difference in the anisotropy of the 21-cm signal. 

In all of our reionization scenarios except one, we 
have assumed that the rate of recombination is uniform 
everywhere in the IGM. However, in reality, the recom¬ 
bination rate is expected to be dependent on the den¬ 
sity of the ionized medium. Specifically, dense structures 
with sizes on the order of a few kpc, mostly unresolv- 
able in this type of simulations, are expected to boost 
the recombination rate significantly and thus increase the 
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Reionization 

scenario 

uv 

UIB 

SXR 

PL 

n 

Non-uniform 

recombination 

Fiducial 

100% 

- 

- 

1.0 

No 

Clumping 

100% 

- 

- 

1.0 

Yes 

UIB dominated 

20% 

80% 

- 

1.0 

No 

SXR dominated 

20% 

- 

80% 

1.0 

No 

UV-tSXR+UIB 

50% 

10% 

40% 

1.0 

No 

PL 2.0 

- 

- 

- 

2.0 

No 

PL 3.0 

- 

- 

- 

3.0 

No 


Table 1. The relative contribution from different source types in 
our reionization scenarios. UV=Ultra-violet; UlB=Uniform ion¬ 
izing background; SXR=Soft X-ray background; PL=Power-law. 


number of ionizing photons required to complete reioniza¬ 
tion. They will also give rise to self-shielded regions like 
Lyman limit systems. There has been some effort in in¬ 
cluding these effects in the simulations of the large scale 
EoR 21-cm signal (Choudhury, Haehnelt & Regan 2009; 
Sobacchi & Mesinger 2014; Choudhury et al. 2015; Shukla 
et al. in prep). We have included the effects of these shelf- 
shielded regions in the IGM by using equation (15) in 
Choudhury, Haehnelt & Regan (2009). Although this ap¬ 
proach somewhat overestimates the impact of non-uniform 
recombinations (due to the coarse resolution of the density 
fields that we use in our semi-numerical models) it can still 
serve to illustrate their effect on the reionization topology. 


2.3 Generating redshift-space brightness 
temperature maps 

Following the steps described in Section 2.1 and 2.2 we have 
simulated coeval volumes of the matter density, ionization 
and velocity fields at various redshifts for the different mod¬ 
els of reionization described in Table 1. The density and ion¬ 
ization fields are then combined to generate the brightness 
temperature maps assuming that Ts 3> Tcmb^, where Ts 
and Tomb are the spin temperature and the CMB temper¬ 
atures respectively. We then take into account the peculiar 
velocities in order to construct the redshift-space signal. By 
redshift space, we mean the space that will be reconstructed 
by an observer assuming that all redshifts are purely due to 
the Hubble expansion. An emitter with a line-of-sight pecu¬ 
liar velocity U|| at a real-space position r will be translated 
to an apparent, redshift-space position s following 


® Note that we have not taken into account the effect of fluc¬ 
tuations in the spin temperature. Spin temperature fluctua¬ 
tions due to Lyman-a pumping and heating by X-ray sources 
can affect the 21-cm brightness temperature fluctuations sig¬ 
nificantly during the early stages of the EoR (Mao et al. 2012; 
Mesinger, Ferrara & Spiegel 2013; Ghara, Choudhury & Datta 
2015). However, as we discuss in the later parts of this paper 
(Section 4.1), the anisotropy in the 21-cm power spectrum will 
not be strongly affected by this, especially once the early phase 
of EoR is over (i.e. when xjj, < 0.95). 


For a more in-depth explanation of 21-cm redshift-space dis¬ 
tortions, see e.g. Mao et al. (2012), Jensen et al. (2013) or 
Majumdar, Bharadwaj & Choudhury (2013). 

We implement the redshift-space distortions using the 
same method as in Jensen et al. (2013). This method splits 
each cell into n smaller sub-cells along the line-of-sight and 
assigns each sub-cell the brightness temperature STt(r)/n, 
where STh(r) is the temperature of the host cell. It then 
moves the sub-cells according to Equation (2), and regrids 
them to the original resolution. Here, we use n — 50, 
which gives results that are accurate down to approximately 
one-fourth of the Nyquist wave number, or k < kt^/4 = 
7r/4 X 600/(714 Mpc) = 0.66 Mpc“^ (Mao et al. 2012; 
Jensen et al. 2013). Figure 2 shows the resulting redshift- 
space brightness temperature maps for six of the seven dif¬ 
ferent reionization scenarios considered here (at a time when 
the reionization had reached 50%). 


3 QUANTIFYING THE POWER SPECTRUM 

ANISOTROPY 

Since only the line-of-sight component of the peculiar ve¬ 
locities affects the 21-cm signal, redshift-space distortions 
will make the observed 21-cm power spectrum anisotropic. 
To quantify this anisotropy, it is common to introduce a 
parameter /r, defined as the cosine of the angle between 
a wave number k and the line-of-sight. The redshift-space 
power spectrum of the EoR 21-cm signal at large scales 
can be written as a fourth-order polynomial in p under 
the quasi-linear approximation for the signal presented by 
Mao et al. (2012). This quasi-linear model is an improve¬ 
ment on the linear model of Bharadwaj & Ali (2005) and 
Barkana & Loeb (2005) for the 21-cm signal. In the quasi- 
linear model one assumes the density and velocity fluctua¬ 
tions to be linear in nature but the fluctuations in the nentral 
fraction are considered to be non-linear. The expression for 
the 21-cm power spectrum, P®(fc,/i), becomes: 

F^(k,/.) = m'(z) [PpH.,PH.(fc) + 

+ 2^ ^PHi.Pm(^) + U ^pm.pm(^)] 1 (3) 

where pni is the neutral hydrogen density, and pM is the 
total hydrogen density. 

It has been suggested that this form of de-composition 
can be used to extract the cosmology from the astro¬ 
physics on large spatial scales, by extracting the coefRecient 
of the p^-term in Equation (3) (Barkana & Loeb 2005; 
McQuinn et al. 2006; Shapiro et al. 2013). However, this is 
very challenging due to the high level of cosmic variance at 
these length scales and also due to the fact that this decom¬ 
position is in a non-orthonormal basis. 

However, the astrophysical information contained in 
the term is also interesting. This term is deter¬ 
mined by the cross-power spectrum of the neutral mat¬ 
ter density and the total matter density, which is a mea¬ 
sure of the “inside-outness” of the reionization topology 
(Majumdar, Bharadwaj & Choudhury 2013; Jensen et al. 
2013; Ghara, Choudhury & Datta 2015). However, extract¬ 
ing this term from noisy data using Equation (3) is also 
challenging since power tends to leak over between the p^ 
and p^ terms (Jensen et al. 2013). 

One way to deal with this is to measure the sum of the 
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Figure 2. Redshift space brightness temperature maps for six of our reionization scenarios. The line of sight is along the y-axis. The 
mass averaged neutral fraction for all the panels is xjji = 0.5. 


fj? and fi'^ terms (Jensen et al. 2013). This quantity is much 
more resistant to noise. It will contain a mixture of cos¬ 
mological and astrophysical information but since the cos¬ 
mological matter power spectrum evolves only slowly and 
monotonically, the evolution of the sum of the p terms will 
be mostly determined by the term. 

A different approach is to instead expand the 
power spectrum in the orthonormal basis of Leg¬ 
endre polynomials—a well known approach in the 
held of galaxy redshift surveys (Hamilton 1992, 1998; 
Cole, Fisher & Weinberg 1995). In this representation, the 
power spectrum can be expressed as a sum of the even mul¬ 
tipoles of Legendre polynomials: 

P^(A:,p)= ^ P(^*)Pr(fc). (4) 

I even 

From an observed or simulated 21-cm power spectrum, one 
can calculate the angular multipoles Pf: 

Ptik) = I P{fi)P^ik)d^- ( 5 ) 

The integral is done over the entire solid angle to take into 
account all possible orientations of the k vector with the line- 
of-sight direction. The estimation of each multipole moment 
through Equation (5) will be independent of the other, as 
this representation is in an orthonormal basis. 

Under the quasi-linear model of Mao et al. (2012), only 
the hrst three even multipole moments will have non¬ 


zero values (Majumdar, Bharadwaj & Choudhury 2013; 
Majumdar et al. 2014): 


Po = sn 


P| = 45T, 


ip _|_ p +-P 

PM.PM V 1 PHI.PHI V gJ PHI.PM 


ip +ip 

PM>PM ^ PHnPM 


Pi — 35 ^^^ Ppyi,PM 


( 6 ) 

( 7 ) 

( 8 ) 


The monopole moment Pq is, by dehnition (Equation 5), 
the spherically averaged power spectrum. We see that the 
quadrupole moment, P|, is a linear combination of Pp^.PM 
and PpHi.PM I ®um of the and terms in Equa¬ 

tion (3). Both the sum of the and terms in Equation 
(3) and PI in Equation (7) thus contain the same physical 
information. However, in contrast to the decomposition of 
the power spectrum in terms of powers of p, the expansion in 
the Legendre polynomials is in an orthonormal basis, which 
means that the uncertainty in the estimates of one moment 
will be uncorrelated with the uncertainty in the estimates 
of the other moments. Therefore, the multipole moments 
should be easier to extract from noisy data. This also implies 
that, if measured with statistical signihcance, each multipole 
moment can be used as an independent and complementary 
estimator of the cosmological 21-cm signal. 

To test this, we took simulated brightness temperature 
volumes from the fiducial model (see Section 2.2) and added 
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Figure 3. Comparison of the two different ways of measuring 
the anisotropy. The top panel shows the anisotropy reconstructed 
from noisy data (see the text for details) as a function of neutral 
fraction. The error bars show the 2a variation for 50 different 
noise realizations. The bottom panel shows a histogram of the 
reconstructed noisy values divided by the true values. A narrower 
distribution around 1 is an indication of a better reconstruction. 


50 different realizations of Gaussian noise (later on, we 
will make predictions with realistic noise for LOFAR 

and the SKA). After calculating the power spectra of the 
noisy data volumes, we fit a fourth-order polynomial in /r 
using standard least-squares fitting, and also estimated the 
P| moment using Equation (5). The results of this test are 
shown in Figure 3. The top panel shows the measured values 
of the anisotropy using the two methods. The solid lines 
show the true, noise-free values, and the error bars show the 
2(7 variation in the noisy measurements across the different 
noise realizations. In the bottom panel, we show a histogram 
of the ratio between the measured and the true value for two 
methods. While the difference is not dramatic, we see here 
that the P| is slightly more resistant to noise than the sum 
of the fj, terms. Because of this and also due to the fact that 
each multipole moment represented in the orthonormal basis 
of of Legendre polynomials will be independent of the other, 
we will use the quadrupole moment Pi, as our measure of 
power spectrum anisotropy for the rest of this paper. 

Fluctuations in spin temperature Ts, introduced due to 
the heating by the very early astrophysical sources will af¬ 
fect the monopole moment (Pq) or the spherically averaged 
power spectrum significantly (Fialkov, Barkana & Visbal 
2014; Ghara, Choudhury & Datta 2015). However we expect 


that they will not have much impact on the anisotropy term 
or the quadrupole moment (PI)- We direct the reader to 
Section 4.1 for further details. 

3.1 How good is the quasi-linear model for 
interpreting the anisotropy? 

The quadrupole moment gives us a measure of the power 
spectrum anisotropy. Ideally, we would want to connect 
this measure to the topology of reionization. The quasi- 
linear approximation gives us a way to do this in the early 
stages of reionization and for large length scales. When this 
approximation holds. Equation (7) tells us that P| is a 
weighted sum of the matter power spectrum and the cross¬ 
power spectrum of the total and neutral matter densities. 
This cross-power spectrum measures how strongly corre¬ 
lated the Hi and matter distributions are. If reionization 
is strongly inside-out, this quantity will decrease quickly as 
reionization progresses, and take on a negative value (see e.g. 
Jensen et al. 2013; Majumdar, Bharadwaj & Choudhury 
2013; Ghara, Choudhury & Datta 2015). 

To some degree, the success of the quasi-linear approx¬ 
imation depends on the source model. In Figure 4, we show 
the ratio between the true value of P| and the quasi-linear 
expectation, constructed by first calculating Pp^.p^ and 
Ppjji^PM and then combining them according to Equation 
(7). We see that at k = 0.12 Mpc“^, the approximation 
works rather well (i.e. the ratio is close to 1) for all scenarios 
up to around xm ~ 0.5. For larger spatial scales, the effect 
of sample variance becomes significant, while at smaller spa¬ 
tial scales, non-linearities start influencing the results. The 
quasi-linear approximation for the signal works poorly for 
the PL 2.0 scenario, but works almost perfectly until very 
late stages of EoR for the UIB Dominated scenario. Around 
xhi = 0.7, the quadrupole moment crosses zero, giving rise 
to sudden jumps in the ratio. 

Note that the quadrupole moment of the power spec¬ 
trum is a quantity that can always be calculated from the 
data, even when the signal is highly non-linear and also re¬ 
gardless of the assumed model for the signal. The quasi- 
linear model is simply a tool to help with interpreting the 
quantity. The quadrupole moment will still be useful as an 
independent estimator of the signal and as an observable 
to distinguish between different source models even in the 
non-linear regime. Here, all we can do is look at how this 
estimator evolves in simulations with different reionization 
scenarios, which we discuss in the next section. 


4 REDSHIFT-SPACE DISTORTIONS IN 

DIFFERENT REIONIZATION SCENARIOS 

In Figure 5, we show the monopole moment of the power 
spectrum, Pq —which is by definition the spherically aver¬ 
aged power spectrum—for the reionization scenarios listed 
in Table 1. From here onwards we show all our results for the 
wave number k — 0.12Mpc“^. It is expected that LOFAR 
will be most sensitive around this length scale (Jensen et al. 
2013) and also the quasi-linear approximations seems to 
work quite well for this length scale (Figure 4). The quantity 
Pq has been studied extensively before (e.g. Santos et al. 
2010; Mesinger, Furlanetto & Gen 2011; Mao et al. 2012; 
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Figure 4. The ratio between the quadrupole moment measured from the simulated 21-cm signal, and the quasi-linear approximation, 
calculated from Equation (7) for our different reionization scenarios. 



Xhi 


Figure 5. The monopole moment of the power spectrum for all 
of our reionization scenarios as a function of the global neutral 
fraction at fc = 0.12 Mpc“^. 


Majumdar, Bharadwaj Choudhury 2013; Jensen et al. 
2013; Ghara, Choudhury & Datta 2015 etc.), and the trends 
we are observing here mostly agree with these previous stud¬ 
ies. 

At the large scales we consider here, the behaviour of 
Pq can be readily understood by the quasi-linear approxi¬ 
mation (Equation 6). Initially, the brightness temperature 
fluctuations mostly come from fluctuations in the H l distri¬ 
bution, and Pq decreases as the density peaks are ionized. 
Later on, at a neutral fraction around 0.3 — 0.5 (depend¬ 
ing on the reionization scenario), a peak in Pq appears, as 
large ionized regions form and become the main contribu¬ 
tors to the brightness temperature fluctuations. This peak 
is especially pronounced in the PL scenarios, where the ion¬ 
ized bubbles are particularly large. However, it is completely 
absent in the UIB Dominated scenario, since due to the pres¬ 
ence of a uniform ionizing background the Hi distribution 
at all stages follows the product of the matter distribution 
and the mean brightness temperature, which causes Pq to 
fall off monotonically with decreasing xhi. 

Figure 6 shows the quadrupole moment, P|i which we 
use as a measure of the power spectrum anisotropy. It was 
also studied in Majumdar, Bharadwaj & Choudhury (2013) 
and Majumdar et al. (2014) for EoR scenarios similar to 
the fiducial and clumping scenarios described in this pa¬ 
per. As we saw previously (Section 3.1), the behaviour of 
P| can be understood to some degree using the quasi-linear 
model (Equation 7), where PI is a linear combination of the 
matter power spectrum Ppj^,p^ and the cross-power spec- 
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Figure 6. The quadrupole moment of the power spectrum for all 
of our reionization scenarios as a function of the global neutral 
fraction at fc = 0.12 Mpc“^. 

trum between the total matter density and the neutral mat¬ 
ter density, Pphi.pm- We show the product of the square 
of the mean brightness temperature {ST^, ) and the cross¬ 
power spectrum for the different reionization scenarios in 
Figure 7, along with the product of the STh and the mat¬ 
ter power spectrum (which is the same for all the scenar¬ 
ios). Since the matter power spectrum increases slowly and 
monotonously and JTb decreases rather rapidly with the de¬ 
creasing xhi, the evolution of the term STh Ppm,pm driven 
by the evolution of STh ■ As the contribution from the term 
STh PpM.PM i® the same for all the reionization scenarios 
considered here, the major differences between P| for differ¬ 
ent reionization scenarios come from the differences in the 
evolution of STh PpM,pui - Going back to Figure 6, for most 
scenarios P| initially increases in strength and then starts 
falling. This is because when the IGM is completely neutral, 
Phi = Pm and PpuuPM '"'ill increase when the matter fluc¬ 
tuations grow. As the most massive peaks are ionized, how¬ 
ever, PH I becomes anti-correlated with pM and P| becomes 
negative. The stronger the anti-correlation—i.e. the more 
inside-out the reionization topology is—the more negative 
P| becomes. The scenarios that stand out from the rest are 
again the UIB Dominated scenario, where the reionization 
is more outside-in than inside-out, and the PL 3 scenario. 
In the latter, reionization is driven mainly by a few massive 
sources that produce very large regions of ionized hydrogen 
(Hll). Since these Hll bubbles are so large, the correlation 
between the total matter and H l field becomes very weak at 
large scales, and so PI remains close to zero for most of the 
reionization history. 

Since P| and P| are two independent but complemen¬ 
tary measurements of the 21-cm signal, one can thus visu¬ 
alize the evolution of the 21-cm signal in a reionization sce¬ 
nario as a trajectory in the phase space of P| and PI. In Fig¬ 



Figure 7. The matter and H i density cross-power spectrum as 
a function of the global neutral fraction for all the reionization 
scenarios at fc = 0.12 Mpc~^. The black solid line shows the 
evolution of the matter power spectrum with xh i > which is same 
for all the reionization scenarios considered here. 

ure 8 we show the reionization scenarios considered here in 
such a phase space diagram. Looking first at the fiducial sce¬ 
nario, it initially moves upwards to the right, with both P| 
and P| increasing as the matter fluctuations grow. At a ion¬ 
ization fraction of a few percent, it starts moving downward 
in an arc, eventually reaching another turn-around point at 
xhi ~ 0.4. The point where the Hi fluctuations are at their 
strongest seems to roughly coincide with the point where P| 
reaches its minimum. After this, the trajectory approaches 
the (0,0) point of the phase space, as the IGM becomes more 
and more ionized and the 21-cm signal starts to disappear. 

From Figures 8 and 6, we see that the power spectrum 
anisotropy (i.e. P|) evolves remarkably similarly in most of 
the reionization scenarios, even when the spherically aver¬ 
aged power spectrum (i.e. P|) differs. The slope is very sim¬ 
ilar in the early stages of reionization, and the minima in 
P| occur at around xhi ~ 0.5 for all reionization scenarios 
except the extreme UIB Dominated scenario. 

In Appendix A we show that the strongest contributing 
factor to the evolution of P| is the phase difference between 
the matter and the H l fields. As long as a major portion 
of the ionizing flux originates from the high density regions 
(or the collapsed objects), the phase difference between the 
two fields will evolve in roughly the same way regardless of 
the strength of the H l fluctuations. Only in extreme scenar¬ 
ios such as our UIB Dominated scenario—where the flux is 
evenly distributed—does the phase difference deviate from 
the other scenarios. 

Thus, the power spectrum anisotropy will probably not 
be useful in telling different source models apart, unless the 
sources are extreme in some sense. However, it appears to 
offer a robust way of measuring the history of reionization. 
We discuss this in more detail in the next section. 
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Figure 8. Trajectories of the different reionization scenarios in 
Pj — Pq phase space for k = 0.i2 Mpc“^. The corresponding 
values of xjj i have been printed on the trajectory of the fiducial 
scenario show the global neutral fraction. For clarity, this has 
been shown only for the fiducial scenario. 


4.1 The effects of spin temperature fluctuations 


In the reionization scenarios discussed here we have not 
considered the effect of spin temperature fluctuations due 
to Lyman-Q pumping and heating by X-ray sources. These 
effects can influence the 21-cm brightness temperature fluc¬ 
tuations significantly during the early stages of EoR, if the 
heating of the IGM is late (Fialkov, Barkana & Visbal 
2014). This can also affect the line-of-sight 
anisotropy in the signal (Ghara, Ghoudhury & Datta 
2015; Fialkov, Barkana & Cohen 2015; 

Ghara, Datta & Ghoudhury 2015). In a late-heating 
scenario, if we define the spin temperature fluctuations in 
the Hi distribution as 


77 ( 2 :, x) = 1 - 


Tcmb{z) 
Ts{z,x) ’ 


( 9 ) 


then the angular multipole moments of the redshift space 
power spectrum under the quasi-linear approximations will 
take the form 


Po = sn (z) 

+ ‘^Pri,PHi 


PM«PM 




PHi.PHi 


+ Fn,TJ + 


2 2 
P g^PHl.PM "I" g^»7,PM 


Pi = A STb (z) 


^PM.PM P ^PpHl,PM P ^PptPM 


Pl = ^&Tb («)PpM,PM 


( 10 ) 

( 11 ) 

( 12 ) 


As shown by Fialkov, Barkana & Visbal 
(2014); Ghara, Ghoudhury & Datta (2015) and 
Ghara, Datta & Ghoudhury (2015), the spin temperature 
fluctuations affect the monopole moment of the power spec¬ 
trum (i.e. the spherically averaged power spectrum) most 


severely. Figures 8 and 9 in Ghara, Ghoudhury & Datta 
(2015) show that the amplitude of the power spectrum gets 
amplified significantly at all length scales (by at least two 
orders of magnitude) during the cosmic dawn (when the 
first X-ray sources emerge) and there is a prominent dip in 
its amplitude at large length scales during the early phases 
of reionization (i.e. when Hhi ~ 0.95). This behaviour of the 
spherically averaged power spectrum is due to the fact that 
during the cosmic dawn and the early stages of the FoR, the 
auto power spectrum of the spin temperature fluctuations, 
P,j,r], makes the dominant contribution to the monopole 
moment (see figure 10 of Ghara, Ghoudhury & Datta 2015). 
This term will be zero when one assumes that Ts Tcmb- 
They have not estimated any higher order multipole 
moments of the power spectrum from their simulations. 

However, as we can see from Equation (11), the Pri^n 
term does not contribute to the quadrupole moment which 
we use in this paper to quantify the redshift space anisotropy 
in the power spectrum. The only contribution from spin 
temperature fluctuations to Pi is the cross-power spec¬ 
trum of the spin temperature fluctuations and the mat¬ 
ter density fluctuations, Pri,pyi- In the case of late heat¬ 
ing, Ghara, Ghoudhury & Datta (2015) find that the Prj,p^ 
term is comparable to Pphi.pm large scales during the 
very early stages of the EoR (i.e. xhi ~ 1-0 — 0.95), 
but later on, it becomes few orders of magnitude smaller. 
The evolution of Pr],pyi with redshift and xhi as shown in 
Ghara, Ghoudhury & Datta (2015), may vary depending on 
the properties of the heating sources. These results thus sug¬ 
gest that the quadrupole moment is not completely immune 
to spin temperature fluctuations, but it is much less sensitive 
than the monopole moment. 


5 OBSERVABILITY OF THE 

REDSHIFT-SPACE ANISOTROPY 

Having seen how the power spectrum anisotropy behaves un¬ 
der idealized conditions, we now explore some of the compli¬ 
cations one will encounter when attempting to observe this 
effect. 


5.1 The lightcone effect 

In the previous section, we showed the 21-cm power spectra 
estimated from so-called “coeval’ data volumes, i.e. the di¬ 
rect outputs from simulations, where the signal is at a fixed 
evolutionary stage in the entire volume. Because of the finite 
travel time of light, an observer will not be able to observe 
the signal like this. Instead, the signal will be seen at a dif¬ 
ferent evolutionary stage at each observed frequency. This 
is known as the lightcone effect. 

Because of the lightcone effect, a three-dimensional 
measurement of the 21-cm signal can never be done at a fixed 
single neutral fraction. To measure, say, the power spectrum 
from a real observation, one has to average over some fre¬ 
quency range (corresponding to a range in neutral fraction). 
This averaging has been shown to have a small effect on the 
spherically averaged power spectrum (Datta et al. 2012)'*. 


* The lightcone averaging could have a significant impact on the 
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Figure 9. The monopole and quadrupole moments of the power spectrum for three different reionization scenarios as a function of 
global neutral fraction at fc = 0.12 Mpc“^ estimated from the light cone cubes with different frequency bandwidths. Note that for the 
light cone cube results the xhi values quoted in the x-axis of these plots are the effective xhi values within the sub-volume limited by 
the bandwidth. The shaded regions in three different shades of grey, dark to light, show the la uncertainty in these estimates due to the 
system noise, measured in a bandwidth of 20 MHz, after 1000 and 3000 hours of LOFAR and 100 hours SKA observations, respectively. 


It has also been shown that the lightcone effect does not 
by itself introduce any significant additional anisotropy in 
the 21-cm signal at measurable scales (Datta et al. 2014). 
However, the averaging will have some effect on the mea¬ 
surements of the anisotropy that is already present due to 
the redshift-space distortions. 

In Figure 9, we show the monopole and the quadrupole 
moments of the power spectrum estimated from lightcone 
volumes generated as described in Section 2.3. We com¬ 
pare the results for coeval data volumes with bandwidths of 
10, 15, 20 and 25 MHz for the fiducial reionization scenario, 
along with the two most extreme scenarios: the UIB Domi¬ 
nated and PL 3 scenarios. In general, a narrower bandwidth 
means that the effective averaging is done over a smaller 
range of neutral fractions, and the value of P| and Pq will 
follow the coeval values more closely. However, a narrow 
bandwidth will increase the sample variance. 


spherically averaged power spectrum, especially during the early 
stages of EoR, if one considers the spin temperature fluctuations 
in the signal due to the inhomogeneous heating by the X-ray 
sources and considers averaging the signal for a significantly large 
frequency bandwidth (Ghara, Datta & Choudhury 2015). 


Figure 9 shows that for the fiducial scenario, the light¬ 
cone effect changes P| only marginally, and only at low neu¬ 
tral fractions. The effect is strongest for the other two sce¬ 
narios, but this is mostly due to sample variance, as the sig¬ 
nal is being measured over a narrow slice of the full volume. 
The lightcone effect on Pq is largest around xhi ~ 0.8 for 
the fiducial model. For the other two reionization scenarios 
it is more prominent at low neutral fractions. 


5.2 Detector noise 

Any radio interferometric measurement of the 21-cm signal 
from the EoR will face several obstacles, including detec¬ 
tor noise, galactic and extragalactic foreground sources and 
ionospheric disturbances. While a full treatment of all these 
effects is beyond the scope of this paper, we do investigate 
the fundamental limitations from detector noise and fore¬ 
grounds for LOFAR and the SKA. 

To calculate the detector noise, we use the same method 
as in Jensen et al. (2013). We begin by calculating the u,v 
coverage for the antenna distribution in the LOFAR core 
(Yatawatta et al. 2013). For the SKA, we assume a Gaus¬ 
sian distribution of antennas within a radius of around 2000 
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Figure 10. Effects of foreground subtraction on P|. The figure 
shows the difference between measured after and before adding 
and subtracting foreground for a few different integration times 
for LOFAR and the SKA. The results are shown at fc = 0.12 
Mpc“^ for the fiducial scenario, with the measurements carried 
out at 15 MHz bandwidth. 


meters (Dewdney et al. 2013), assuming 150 antenna sta¬ 
tions. 

We fill the u, v planes with randomly generated Gaus¬ 
sian noise with a magnitude calculated using the formalism 
from McQuinn et al. (2006). We use the same parameters 
for LOFAR as in Jensen et al. (2013). For the SKA we as¬ 
sume a total collecting area of 5 x 10® m^. This collecting 
area, along with the assumptions for the baseline distribu¬ 
tion, is consistent with the current plans for the first-phase 
version of the low-frequency part of the SKA (SKAl-LOW) 
(Koopmans et al. 2015). 

Finally, we Fourier transform the noise in the u, v plane 
to get the noise in the image plane. We do this for many 
different frequencies to get a noise lightcone with the same 
dimensions as our signal lightcones. Here, we do not take 
into account the frequency dependence of the u, v coverage. 

Using these simulated noise lightcones, we can then cal¬ 
culate the power spectrum error due to detector noise. In 
Figure 9, we show these errors for 1000 and 3000 hours of 
LOFAR observations and 100 hours of SKA observations as 
shaded regions (in three different shades of grey, dark to 
light, respectively). We see that for LOFAR, the noise error 
for P| is much higher than the uncertainty due to the light- 
cone effect. Clearly, LOFAR will require a long integration 
time to measure the evolution in P|. For the SKA, however, 
the noise error is almost negligible, even for relatively short 
integration times. In this case, the lightcone effect also be¬ 
comes important, at least at the later stages of reionization. 


5.3 Foregrounds and PSF effects 

In addition to detector noise, there are many other issues 
that will complicate the extraction of P| from a real mea¬ 
surement. One major obstacle is foreground emission from 
Galactic and extra-Galactic sources. To properly measure 
the power spectrum anisotropy, the foregrounds need to be 
removed from the data using some algorithm. 

We simulate foreground emission from diffuse and local¬ 
ized Galactic synchrotron radiation. Galactic free-free radi¬ 
ation and extragalactic bright radio sources using the meth¬ 


ods described in Jelic et al. (2008, 2010). After adding these 
sources to the simulated signal, we use the GMCA algo¬ 
rithm to model and remove them (see Chapman et al. 2013 
for a detailed description). GMCA (Generalized Morpholog¬ 
ical Component Analysis) works by attempting to express 
the full signal+foregrounds in a wavelet basis where only a 
few components are non-zero. Since the signal is very weak 
compared to the foregrounds, it will tend not to be included 
in the wavelet basis and can thus be separated from the 
foregrounds. 

A second complication is the beam, or the point-spread 
function (PSF), of the instrument. The effect of the PSF 
is essentially a smoothing in the sky plane, de-emphasizing 
large k± modes and introducing an anisotropy in the sig¬ 
nal. To properly study the feasibility of extracting P| from 
real measurements, we would need to model the frequency- 
dependent PSF and compensate for its effects, which is be¬ 
yond the scope of this paper. Instead, we simply calculate 
the change in PKJIhi) after adding and subtracting fore¬ 
grounds. This change is shown in Figure 10 for a few dif¬ 
ferent integration times. Here, P^’^® denotes the quadrupole 
moment after adding and subtracting foregrounds, while PI 
is the quadrupole with no foreground effects included. We 
show only the absolute error since the relative error fluctu¬ 
ates wildly where the signal crosses zero. 

We see that for low noise levels, the foregrounds change 
the measured P| by a few mK^. As seen in Figure 6, the 
signal itself ranges between approximately —5 and 10 mK^. 
For the SKA, the foreground removal is almost perfect at 
high values of Shi, while LOFAR observations appear to re¬ 
quire long integration times to properly measure P|. How¬ 
ever, the GMCA algorithm used here is not the only method 
for removing foregrounds, and it remains possible that other 
algorithms perform better for this particular task. 


5.4 Reconstructing the reionization history 

As we saw in Figure 6, the quadrupole moment evolves 
very predictably as a function of xm; much more so than 
the monopole moment or the spherically averaged power 
spectrum (Figure 5). The only scenarios that deviate from 
the rest are the UIB dominated scenario and—for low neu¬ 
tral fractions—the PL 3.0 scenario. A real-life measurement 
would give P| as a function of z, so if we assume that PK^hi) 
is known, we can use this to reconstruct the reionization his¬ 
tory, i.e. Here, we attempt to do this for our mock 

observations. 

We begin by parametrizing xhi{z) as follows: 


XHi{z) 


0 ii Z < Zr 

1 — exp[—/I(z — Zr)] ii Z > Zr 


(13) 


The parameters Zr and /3 give the endpoint and extent of 
reionization respectively (higher values of /3 indicate a more 
rapid reionization). Given values of Zr and /3 and a mea¬ 
surement of Pi{z), we can use Equation (13) to construct 
PI(®Hi). Since we assume that this latter quantity is known, 
we may tune Zr and /3 until the measured P|(a;Hi) matches 
the true one. 

If a measurement of P 2 ’™'’‘‘'’( 2 ) hris Gaussian noise with 
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Figure 11. Fitting the reionization history to match the slope of P| as a function of xni. Left column: Likelihood contours for 1000 and 
3000 hours integration time with LOFAR and 100 hours of integration time with SKA, for different reionization scenarios. The contours 
are located at the value where the likelihood is a factor lower than the maximum value, corresponding to roughly 90 per cent of the 
full distribution. The crosses show the best-fit values of the parameters Zr and /3 (in the absence of noise), and the red circles show the 
true values. Middle column: The reionization histories corresponding to the best-fit values from the left column, compared to the true 
histories. Right column: Same as the left column, but assuming that Zr is known. The dotted vertical lines show the true values of /3. 


an amplitude of a, the likelihood of measuring this value is: 


TT 1 

/ pSjineas r)S,true\2" 

“ [^2,i ~ ^2,i ) 

1 ^ 

2cr2 


where the product is calculated over all redshifts where there 
are measurements available. We are thus looking for the val¬ 
ues of Zr and that give rise to the xhi history that maxi¬ 
mizes the likelihood. 

Judging by Figure 6, P|(:ehi) can be assumed to be rea¬ 
sonably well-known, especially in the early stages of reion¬ 
ization. Here, we therefore only use the part of PI where 
xhi > 0.75, since the slope of P|(xm) is very consistent from 
scenario to scenario. We take the values from the fiducial 
model lightcone to be the “true” P|(a;Hi). More specifically, 
we use a linear fit to the values of PKxhi) for xhi > 0.75 
plotted in Figure 6 to be our Pj 


Figure 11 shows the results of fitting the reionization 
history, jg here taken from our simulated reioniza¬ 

tion scenarios (see Figure 6), and the noise a is calculated as 
described in the previous section. In the left panel we show 
the likelihood contours for LOFAR noise for 1000 and 3000 
hours of integration time and for SKA noise for 100 hours. 
The crosses show the parameter values that maximize the 
likelihood, while the red circles show the true values, i.e. a fit 
of Equation (13) to the true reionization history (since our 
simulations only go down to Hhi ~ 0.2, the value of Zr rep¬ 
resents an extrapolation; c.f. Figure 1). The error contours 
do not explicitly include the effects of sample variance®. The 


® The inherent non-Gaussianity in the EoR 21-cm signal can 
make the cosmic variance at these length scales deviate from 
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middle panel of Figure 11 shows the reconstructed reioniza¬ 
tion histories corresponding to the best-fit parameters in the 
left panel. 

The fact that the best-fit values are very close to the 
true values for all reionization scenarios shows that the evo¬ 
lution of P| can indeed be used as a robust measurement of 
the reionization history. For LOFAR, it is clear that noise is 
an important obstacle. With 1000 hours of integration time, 
we can only put rather loose constraints on the reionization 
history. After 3000 hours, however, the constraints become 
much better. For the SKA, the noise is so low that even for 
100 hours integration time, the reionization history can be 
determined almost perfectly. 

We also see that for our parametrization of the reioniza¬ 
tion history, there is quite a bit of degeneracy between the 
two parameters, as seen by the elongated likelihood con¬ 
tours. In a real-world scenario, it is likely that there will be 
some information available on one of the parameters. For ex¬ 
ample, Zr may be estimated by looking for the redshift when 
the spherically averaged power spectrum tends to zero. In 
the right column of Figure 11, we show the results of fit¬ 
ting only /3, assuming that Zr is known. In this case, LO¬ 
FAR measurements may be able to say something about the 
reionization history, or at least exclude meaningful regions 
of parameter space. 

While all of our reionization scenarios are constructed 
to have very different topologies, they are also tuned to have 
the same reionization history. A widely accepted view is that 
the 21-cm power spectrum is much more sensitive to xhi 
than z (e.g. Iliev et al. 2012), and so we do not expect the 
reionization history to have a large effect on our results. 
Nevertheless, to test whether a different reionization history 
would change the shape of P|(a;Hi), we also constructed a 
scenario with the same source properties as the fiducial sce¬ 
nario, but with the source efhciencies increeised, in order to 
produce an earlier reionization. The results for fitting the 
reionization history for this scenario is shown in the bot¬ 
tom row of Figure 11, still using P|(a;Hi) from the fiducial 
scenario as the true value. While a full investigation of the 
effects of the exact reionization history is beyond the scope 
of this paper, it is clear from this test that a modest shift in 
the history has no significant effect on the results. 

The analysis presented in this section is a proof-of- 
concept of the possibility of extracting the reionization his¬ 
tory from the redshift-space anisotropy of the 21-cm signal. 
In actual radio-interferometric observations, the presence of 
residuals from foregrounds in the reduced 21-cm data can 
introduce further uncertainties in these estimations, which 
we have not accounted for. 


6 SUMMARY AND CONCLUSIONS 

The velocities of matter distort the 21-cm signal and cause 
the power spectrum to become anisotropic. The anisotropy 
depends on the cross-correlation between the density and the 
Hi fields and therefore carries information both about the 


its generally assumed Gaussian behaviour (Mondal et al. 2015; 
Mondal, Bharadwaj &; Majumdar 2015), increasing the uncer¬ 
tainties in the reconstruction of the EoR history. 


cosmology and about the properties of the sources of reion¬ 
ization. In this paper we have explored a number of different 
reionization scenarios with different types of sources of ion¬ 
izing photons in order to study the effect of the sources of 
reionization on the power spectrum anisotropy. 

We hnd that the power spectrum anisotropy, as mea¬ 
sured by the quadrupole moment of the power spectrum, 
P|, as a function of global neutral fraction evolves in a very 
similar way for all reionization scenarios in which the sources 
and the matter distribution are strongly correlated. Only the 
most extreme scenarios, such as reionization driven by a uni¬ 
form ionizing background or by very massive sources, show 
significant deviations from the others. Thus, PI is not a sen¬ 
sitive tool for separating reionization scenarios. However, it 
can be used as a robust tracer of the reionization history. 

Interferometers such as LOFAR or the SKA can only 
measure the fluctuations of the 21-cm signal, not the actual 
magnitude of the signal. Therefore, extracting the reioniza¬ 
tion history is not trivial, and can only be done by fitting 
a model to the observed data. Previous studies have sug¬ 
gested that looking for the peak in the spherically-averaged 
power spectrum (e.g. McQuinn et al. 2007; Lidz et al. 2008; 
Barkana 2009; Iliev et al. 2012) or the variance (Patil et al. 
2014) will give the mid-point of reionization. However, 
as seen in this paper and elsewhere (e.g. Mellema et al. 
2006; Iliev et al. 2006, 2007; Watkinson et al. 2015), the 
spherically-averaged power spectrum can in fact vary sub¬ 
stantially depending on the properties of the sources of 
reionization. The quadrupole moment is much more model- 
independent, as long as the reionization topology is inside- 
out. This is likely to be the case, considering that the lat¬ 
est Planck results on the optical depth to Thomson scat¬ 
tering combined with observations of high-redshift galaxies 
strongly favour galaxies as the main driver of reionization 
(Bouwens et al. 2015). 

We have shown that measuring the reionization history 
using the quadrupole moment works well for all reionization 
scenarios considered in this paper, except for the most ex¬ 
treme ones. For first-generation telescopes such as LOFAR, 
these types of measurements will be challenging due to the 
uncertainties imposed by the noise and foreground emission. 
However, it should be possible to use the power spectrum 
anisotropy to exclude certain reionization histories. In any 
case the anisotropy can be used to verify that the observed 
signal is indeed the 21-cm emission from the EoR. For the 
SKA, neither the noise uncertainty nor foregrounds will be 
an issue, and it will be possible to accurately measure the 
reionization history even with a relatively short integration 
time. 

Of course, our treatment of observational effects here 
is simplified and does not take into account complications 
from, for example, the frequency-dependence of the PSF, 
ionospheric effects or calibration errors. Nevertheless, it 
shows that the 21-cm power spectrum anisotropy due to the 
redshift-space distortions has the potential to be a power¬ 
ful and nearly model-independent probe of the reionization 
history. 
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APPENDIX A: BEHAVIOUR OF THE 
QUADRUPOLE MOMENT 

The similarity in the behaviour of the quadrupole moment, 
P 2 (Figure 6), estimated for different reionization scenarios 
may appear unintuitive. However, it can be explained using 
the quasi-linear model (Equation 7), where the major con¬ 
tribution in the evolution of Pi has been ascribed to the 
evolution of Pp^uPM (Figure 7). 

The cross-power spectrum between two fields, defined 
in Fourier space as A = Ae*® and Y = will be 

PxY = Re(A*y) = IA||B| cos((/> — 0). Thus the cross-power 
spectrum between the matter and the Hi helds (Pphi,pm)i 
essentially a product of three quantities: the amplitudes of 
the matter density fluctuations and the Hi density fluctu¬ 
ations and the cosine of the phase difference between these 
two fields. Since the matter density fluctuations are the same 
for all our reionization scenarios, this means that the differ¬ 
ences in PpHnPM be solely due to the differences in the 
amplitude of their Hi density fluctuations and its relative 
phase with the matter density field. 

The top panel of Figure Al shows the evolution of 
the cosine of this phase difference for different reionization 
scenarios estimated through the cross-correlation coefficient 

’'ph,,Pm(^) = ^’pHnPM(*)/V-PpH,.PH,(fc)^'pM.PM(fc)- ^U obvi- 
ous observation from Figure Al, is that the evolution of 
rpHi.pM is very similar in all of the reionization scenarios, 
except for the UIB Dom scenario. This is because as long 
as the spatial distribution of the major reionization sources 
follow the underlying matter distribution, the relative phase 
between the Hi and the matter density fields remains ap¬ 
proximately the same and also evolves similarly. In other 
words, PpHi.PM evolves similarly as long the reionization is 
inside-out in nature. 

The bottom panel of Figure Al shows the evolution of 
the HI density power spectrum fpH,,PHi > which measures the 
amplitude of fluctuations in the Hi field. The square root 
of BpHi.PHi i® the other quantity effectively contributing to 
the evolution of Bphi.pm- tt is evident from this figure that 
the amplitude and the location of the peak in Tph,,phi i® 
significantly different for different reionization scenarios as 
it is more susceptible to the H l topology. In the cross-power 
spectrum, however, it gets tuned by the phase difference 




Thi 


Figure Al. The top and the bottom panels show the cross¬ 
correlation coefficient (or the cosine of the relative phase) be¬ 
tween the matter and the H i density fields, and the H i density 
fluctuation power spectrum FpHi.PHu respectively, for all of our 
reionization scenarios, as a function of the global neutral fraction 
at fc = 0.12 Mpc~^. 


(^phi.pm) term. This is what makes the contribution from 
PpHnPM iri Pi similar in most of the reionization scenarios. 

According to the quasi-linear model (Equation 6), Pq 
also contains a contribution from PpHi.Pivu ^rit in addition to 
that it also contains contribution from the power spectrum 
of the Hi field, PpHi.pnn which as described before, depends 
only on the amplitude of the Hi fluctuations and not on 
its relative phase with density field. The contribution from 
PpHi.PHi makes the monopole moment Pq more sensitive to 
the Hi topology (or the Hll bubble size distribution) and 
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thus more dependent on the specific reionization scenario 
(a visual comparison of Figure 6 with the bottom panel of 
Figure A1 further strengthens this argument). Because of 
this the P| is a more robust measure of the reionization 
history than Pq . 
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